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Abstract 

We propose a quantitative estimate of the charge gap that opens in the one- 
dimensional dimerized Hubbard model at quarter-filling due to dimerization, 
which makes the system efi^ectively half-filled, and to repulsion, which induces 
umklapp scattering processes. Our estimate is expected to be valid for any 
value of the repulsion and of the parameter describing the dimerization. It is 
based on analytical results obtained in various limits (weak coupling, strong 
coupling, large dimerization) and on numerical results obtained by exact diag- 
onalization of small clusters. We consider two models of dimerization: alter- 
nating hopping integrals and alternating on-site energies. The former should 
be appropriate for the Bechgaard salts, the latter for compounds where the 
stacks are made of alternating TMTSF and TMTTF molecules. 
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I. INTRODUCTION 



A large number of compounds exhibit one-dimensional electronic properties that can be 
described by a quarter-filled Hubbard model with some kind of dimerization. For instance, 
the Bechgaard salts {TMTSF)2X and {TMTTF)2X , where X denotes an anion like C104~ , 
PFq~ , Br~ , etc... can be regarded concerning their electronic structure as being essentially 
one-dimensional systems above a crossover temperature of the order of 3OK0. From 
stoichiometry it is known that there are 3 electrons in the HOMO for each pair {TMTSF)2, 
so that the system is 3/4~filled in terms of electrons or quarter-filled in terms of holes. In the 
following we will always use hole notation and consider quarter-filled systems. A reasonable 
description of these properties should be provided by the dimerized Hubbard model defined 
by the Hamiltonian 

H = Ho + H,^t, (1.1) 
where the kinetic part of the Hamiltonian is (see Fig. |l|) 

^0 = -h E i4,.C,+,,^ + h-C.) - t2 E iiaC^+l,. + h-C) . (1.2) 

ieven,(T jodd,(T 

The operators create particles in the HOMO of TMTSF or TMTTF with spin a. We 
have included two hopping integrals ti and t2 (ti > t2) to describe the dimerization along 
the stac kHi. The dispersion of this model is given by e{k) = ^tl + tl + 2tit2 cos k and is 
depicted in Fig. ^. The important parameters are the dimerization gap = 2(ti —12) that 
opens at the Brillouin zone boundary and the total bandwidth W = 2{ti + ^2)- 

To the hopping part of the Hamiltonian we have to add an interaction part to describe the 
correlation between the electrons. For simplicity, we have chosen the form of the standard 
Hubbard model 

H;^t = U^n^nii , (1.3) 

i 

where U is the on-site Coulomb repulsion and rij^o- = cl^^c^^^. 
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This model should also provide a good description of salts like (^^4)2^-^6, where FA 
stands for Fluoranthenyl, which undergo a Peierls transition leading to dimerization along 
the stacks at a temperature of order 190K&i. There is another class of materials however 
for which a modified version of the previous model is more appropriate. These materi- 
als are related to the Bechgaard salts, but the stacks now consist of alternating TMTTF 
and TMTSF molecule^. A minimal model in this case would be a Hubbard model with 
alternating atomic on-site energies, with a kinetic part of the form 

^0 = -iI](cl,aCi+l,a + h.C.) -eo {4,aCi,a ' 4+l,aC^+l,a) , (1-4) 

i,(T ieven,a 

where 2eo is the energy splitting between the HOMO of the TMTSF and TMTTF. The 



dispersion relation is given by e{k) = JeQ + 2t'^ + 2t'^ cos k. It has the same form as in the 



previous case. The dimerization gap is now given by A/) = y^o + and the bandwidth 
hy W = 2^0 • However, although the dispersion relations are the same, the two models 
are different because the interaction part will not be the same when expressed in terms of 
the operators that diagonalize the kinetic part of the Hamiltonian. Let us also note that a 
similar model has been recently proposed by Sudb0 et a/.i as a one dimensional analog of the 
copper oxide layers in the high temperature superconductors. A model where dimerization 
is induced via alternating on-site repulsions has also been studied!. 

It is clear from Fig. ^ that, due to dimerization, the system is effectively half-filled. In 
this case, it is known from the general theory of one-dimensional mo dell'0 with umklapp 
scattering that an on-site repulsion of arbitrary size will open a gap in the charge sector. 
This has been shown explicitly in the case of the Hubbard model for which exact results 
are available from the Bethe ansatz solutionlli. In the general case, where no exact solution 
is known, a quantitative estimate of the gap is not available so far. This is unfortunate 
because such an estimate is necessary to interprete the activated behaviour of the resistivity 
observed at relatively low temperatures in several compounds, which in turn provides an 
estimate of the magnitude of the Coulomb repulsions. An analysis of that sort has already 
been performed for the Bechgaard salts of the TMTTF family on the basis of a preliminary 



determination of the charge gap in the model with alternating hopping amplitudescj. 

In this paper, we propose a quantitative estimate of the gap in the case of the dimerized 
models described above that should be valid in the whole range of parameters. This estimate 
is based on approximations that provide analytical expressions in various limits, and on 
numerical calculations using Lanczos diagonalization of small clusters. As we think that 
the results might be useful to readers that do not want to go through the details of the 
calculation, we start in Sec. II with a detailed account of the results, including analytical 
expressions, tables and curves, that can be used independently of the rest of the paper. Then, 
in Sec. Ill we explain how and for which values of the parameters numerical estimates can 
be obtained on the basis of exact diagonalization of small clusters. The following sections 
are devoted to the approximate expressions that can be obtained analytically in various 
limits: Sec. IV summarizes exact results in a few trivial limits. Sec. V deals with the large 
dimerization case, where the dimerization gap is larger than both the repulsion and the Fermi 
velocity, and Sees. VI and VII with the case of weak and strong coupling, respectively. The 
summary can be found in Section II. Throughout the paper, subsections A will deal with 
the model with alternating hopping amplitudes, and subsections B with the model with 
alternating on-site energies 

II. SUMMARY OF RESULTS 

In Fig. we have illustrated the different regimes in which analytical results have been 
obtained as a function of the model parameters. In this section, we give a summary of these 
results for each model. 

A. Alternating hopping amplitudes 



4 



1. Numerical estimates: 

Using exact diagonalization of small clusters, we have obtained estimates of the charge 
gap for intermediate values of repulsion U and for ^2 ^ 0.5ti. The limitation comes from the 
size of the clusters (maximum 16 sites), so that only large gaps corresponding to correlation 
lengths smaller that 16 sites could be extracted. The numerical results are shown in Fig. ^a. 
They are in good agreement with analytical approximations in different limits, as we will 
show in the subsequent paragraphs. 



2. Large dimerization: 



In this limit we have mapped our model onto the exactly solvable Hubbard model at 
half-filling and we found that the gap is given by: 



^ ^ /f^g-27rt2/C/g-7rt2/4ti 

2 7r\ U 



U W 
I 2 



16ti 



2h 



if ti > t2 > U; 



if ti > f/ > t2. 



(2.1) 



In Fig. 1^ we have compared our estimate with numerical results (curve a) for t2 = 0.1, where 
we can see the crossover from the exponential regime (curve d) to the linear one (curves b 
and c, b is without the quadratic term U^/16ti) for U/ti ~ 0.5 . 



3. Weak coupling limit 

In the weak coupling limit {U <^ ^1,^2), using the RG group method and results from 
the large dimerization limit, we found that the gap is exponentially small [see Eq. ( |6.24| )]: 




where the parameters a and b can be found in Table | for different values of ^2/^1- 
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In Fig. P we compare the numerical values of the gap for t2 = 0.3 and t2 = 0.4 with the 
weak coupling approximation Eq. ( |2.2| ). We can see that up to = 4 the weak coupling 
formula gives very good results. 



4- Strong coupling limit 

When the on-site repulsion is large enough {U ^ ^1,^2? 4t^/t2), the effective Hamiltonian 
of our model is a t — J like Hamiltonian. Applying degenerate perturbation theory, we found 
that the charge gap is equal to the dimerization gap with a 1/U correction [see Eq. (|7.20| )]: 



A^(l-c^) , (2.3) 



U 

where the parameters A^/ti and c can be found for different values of ^2/^1 in Table |[ 

5. Intermediate region 

On the basis of the previous results, we can propose an estimate of the gap for any value 
of the parameters. This can be achieved by using Fade approximants to connect the exact 
results we have obtained for U small and U large respectively. For t2 > 0.2, we have used the 
weak coupling expression up to Uq, where Uq is not too large (actually we choose Uo/ti = b 
for ^2/^1 ^ 0.5 and Uo/ti = 4 for ^2/^1 > 0.5) and Fade approximants of the form 

62 + hiU + 

for U > Uq. The coefficients are chosen in such a way that, for U +00, the Fade 
approximant gives the correct large U behavior given by Eq. ( p.3| ) and that, aX U = Uq, the 
resulting curve and its derivative be continuous. The curves depend weakly on the value 
of Uq\ Changing the value of Uq by 50% affects the value of the gap by less then 10%. In 
Table |II| we give the values of the coefficients and of Uq for different values of 12/ ti. 

For t2 = O.lti we have also used a Fade approximant, but now the limiting behaviors 
were determined by the second equation in ( |2.1|) , which is valid for the intermediate values 
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of U for t2 small, and by Eq. ( |2.3| ) in the large U limit. At the crossing point of the weak- 
coupling expression and the Pade approximant, we have got rid of the kink (see Fig. |^) by 
a linear interpolation. The resulting curves for t2/ti ranging from 0.1 to 0.9 are presented 
in Fig. |. 



B. Alternating on site energies 



1. Large dimerization: 



In this limit we have mapped our model to the exactly solvable Hubbard model at half- 
filling and we found that the gap is given by: 



'8 U 

-t\ — exp{-7rt'^/Ueo) , if eo > t > f/; 

TT V 2^0 



(2.5) 



2f^ 
U- — 



if eo > f/ > t. 



2. Weak coupling limit 



In the weak coupling limit {U ^ ^o,^); on the basis of renormalization group analysis 
and results from the large dimerization limit, the gap is exponentially small and is given by 
[see also Eq. ( |6.38| )]: 



Ac = atJjexp fy-bjj 



(2.6) 



The parameters a and h are given in Table III for some values of eo/t. 



3. Strong coupling limit 

When the on-site repulsion is large enough {U ^ eo, t, 4t^/eo), the effective Hamiltonian 
of our model is again a.t — J like Hamiltonian and we found that the charge gap is [see also 
Eq iVM]: 
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Ac = (l - c^) . (2.7) 
The parameters A^, and c are given in Table |T| for some values of eo/t. 

III. NUMERICAL SIMULATIONS 

We have done intensive numerical simulations based on Lanczos diagonalization of small 
clusters for the model with alternating hopping integrals. Our numerical estimation for the 
gap was based on the usual formula 

Ae = ^ lim A,{N) , 

iV— ++00 

A^{N) = E{N + 1; 2A^) + E{N - 1; 2N) - 2E{N; 2N) , (3.1) 

where E{M] L) is the ground state energy for M particles on L sites. We have obtained 
results for systems of 4, 8, 12 and 16 sites. Periodic (resp. antiperiodic) boundary conditions 
have been used for L=8 and 16 (resp. L=4 and 12) to have open-shell systems. Typical 
results are shown on Fig. ^, where we have plotted Ac{N) as a function of 1/A^. Good 
estimates can be obtained only when the exponential regime Ac(A^) — ~ exp(— cst/A^) 
has been reached for N = 16. But the exponential regime is obtained when Vc/N < Ac, 
that is > Vc/Ac- So we can only obtain good numerical estimates of the gap when Ac is 
not too small, i.e. when U and A^ are not too small. The gap can then be extracted by use 
of the Shanks transformationlii given by 

A(iVi)A(iV3)-A(iV,)^ 
" A{Ni) + A{Ns)-2A{N2) ' ^ ' ^ 

This procedure has been shown to work very well in the case of the Haldane gap of spin-1 

chainJii. Applying it for (Ni, N2, N^) = (4,8,12) and (8,12,16) respectively gave similar 

results. The results quoted throughout this paper are those obtained for (8, 12, 16) because 

they are a priori better. 

Another limitation comes from the size of the repulsion. When U is very large, the 

convergence of the Lanczos algorithm becomes very slow, and it is no longer possible to get 
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good values of the energy. So, in summary, good estimates have been obtained for A/5 ^ ti 
and 1 ^ U/ti ^ 60. This does not cover the range of parameters for actual compounds, 
for which one often has Ao/ti, or U/ti <^ 1, and analytical methods are clearly needed to 
complement these numerical results. In fact, the analytical methods developed in the rest of 
the paper provide an estimate of the gap for any value of the parameters, and the numerical 
simulations have been only used to check the analytical results. 



A schematic picture of the regions of interest as a function of the model parameters 
are given in Fig. ^ for each model. Exact, although trivial, results can be deduced on the 
boundaries and they are summarized in this section. 



(i) ti = 12'. This is a quarter-filled Hubbard model and we know from the Bethe ansatz 
solution that there is no gap in the spectrum. 

(ii) t2 = 0: The Hamiltonian describes a set of independent systems, each system con- 
sisting of two sites. At 1/4 filling there is one electron for each pair of sites in the ground 
state. The gap is then given by Ac = Eq{0) + Eo{2) — 2_E'o(l)! where Eq^ti) is the ground 
state energy of n electron on a pair of sites. These energies are Eq{0) = 0, Eq{1) = —2ti 
and Eo{2) = {U — \Ju'^ + 16ti)/2 and the charge gap is 



IV. EXACT RESULTS 



A. Alternating hopping amplitudes 



r u 



if h > U; 




(4.1) 



(iii) U = 0: The band structure is given on Fig. 0. The Fermi energy is in the middle of 



the lower band, and = 0. 
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(iv) U = +00: The energies that enter Eq. (|3.1|) are the same as for free spinless fermions 
at half-fiUing, because the energy is independent of the spin when U = +00 for open 
boundary conditions (or for periodic one in the hmit +00). So, the charge gap is 

actually the same as for a system of spinless fermions described by the kinetic part of the 
Hamiltonian at half-filling, which is nothing but the dimerization gap (see Fig. So 

= Ad. 

B. Alternating on site energies 

(i) Eq = 0: This is a quarter filled Hubbard model, and there is no gap. 

(ii) Eq = +00: The odd sites are outside the Hilbert space. The ground state has one 
particle per even site, and the first excited state has a doubly occupied site. So clearly 
A, = f/. 

(iii) U = 0: The Fermi energy is in the middle of the band, and Ac = 0. 

(iv) U = +00: The dimerization gap 2eq and the charge gap coincides, like in the case 
of alternating hoppings. 

(v) In the atomic limit {t = 0) the charge gap is given by A^ = min(f/, 2eq). 

V. THE LIMIT OF LARGE DIMERIZATION 

In this section, we consider the limit of large dimerization, that is the limit where the 
dimerization gap A/) is much larger than both the width of each subband and the repulsion 
U. This limit corresponds to ti ^ t2,U for the model with alternating hoppings and to 
£0 ^ t,U for the model with alternating on-site energies. It is particularly interesting 
because in both cases the model can be mapped onto the half-filled Hubbard model, so that 
we can use the exact result provided by the Bethe Ansatz for the charge gapE. 
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A. Alternating hopping amplitudes 



To diagonalize the part of the Hamiltonian corresponding to the hopping terms ti, we 
introduce bonding and antibonding operators defined by: 



1 

'''''' ~ V2 



(5.1) 



where j is now even. In terms of these operators, the Hamiltonian becomes 



J even, a 

t 



2 . 

jeven,cr 

Jeven 

The bonding and the antibonding bands are separated by a large energy 2ti. The occupation 
of the upper band is thus negligible and to zeroth order in 1/ti we get 

^eff = -| E (^L^j+2,a + h.c.) + E ^b,j^nbj,l + 0{l/ti) , (5.3) 

J even, (T J even 

which is nothing but the regular Hubbard model with a repulsion U — U/2 and a hopping 
integral i — t2/2. Here the Hamiltonin acts on a Hilbert space where the antibonding 
states arc all empty. The value of the gap in the half-filled Hubbard model is known 
exactly from the Bethe Ansatz solution. For small interaction {U -C t), it is given by 
Ac = {8t/7r)^/W/iex.p{-27rU/t), which, with our notations, reads 

Ae=|-#e--W.. (5.4) 

For large interaction {U :$> t2 but still U <^ti), the exact expression becomes Ac = C//2— 2^2- 
To compare with other hmits, it is useful to go to next order in l/h. Using a Schrieffer- 
Wolff transformation, we can determine the 1/ti corrections to the zeroth order effective 
Hamiltonian by including scattering processes to the antibonding band and we get 
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t fuu'^\ 

jeven,a \ 1 / j'even 

+ ^ E (^WW + &Ww) + o(iA?)- (5.5) 

The on-site repulsion is reduced by a f/^/16ti, and a second nearest neighbor hopping 

appears. The formula for the gap is now modified to 

/Ee-2-*2/c/e-^*2/4ti -^^^ 
2 vr V ^2 

A. = <! (5.6) 

U f/2 

2t2 , if [/ > t2. 



I 2 16ti 

The second formula agrees with Eq. ( [4.1|) when ^2 = 0. The first one will be used in Section 
VI. We have illustrated the above estimates of the gap on Fig. ^. 

B. Alternating on site energies 

In the limit Eq ^ U,t the occupation of the energetically lower lying even sites is much 
larger than that of the odd sites. Using this, we can again find an effective Hamiltonian 
starting from t = 0, in which case only the even sites are occupied. Switching on the hopping, 
the electrons can hop to the energetically unfavorable odd sites, and from those sites they 
can hop further. This second order virtual process produces an effective hopping of order 
t'^/2eQ between the even sites, so that the effective Hamiltonian is given by 

Hes = -— E (S>S+2,<x + h.c.) + t/ E ^i,T%,i • (5-7) 

jeven, a jeven 

Here the Hamiltonian acts on a Hilbert space with empty odd sites. This effective Hamil- 
tonian describes a half-filled Hubbard model with a hopping amplitude t = /2eQ and 
repulsion U = U, so again we can use the expression of the charge gap for the Hubbard 
model. In the weak couphng limit {U <^ f^/so), we get 



In the intermediate coupling limit, where U is larger than the effective hopping, but smaller 
then the energy splitting of the two bands eo, the gap is given by Ac = f/ — 2t^ /eq. 
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VI. WEAK COUPLING LIMIT 



In this limit we can first diagonalize the hopping part of the Hamiltonian and treat 
the interaction as a perturbation. To do that, we introduce the Fourier transforms of the 
electron creation and annihilation operators keeping in mind that there are two sites in the 
unit cell: 

' J even 

c.,. = ({y.^-''''-'''\.- (6-1) 

' J odd 

Here L is the number of sites, L/2 is the number of unit cells, and the length of the unit 
cell is set to 1. 



A. Alternating hopping amplitudes 



The Fourier transform of the kinetic part of the Hamiltonian ( |1.2| ) is 



(6.2) 



Let us introduce annihilation operators for the electrons in the lower {dk,a) and upper 
(/fc,^) bands by 

1 



di.. 



V2 
1 



where s{k) = expi{ak/2 + k/4) and au is defined by 



(6.3) 



h — ti k 
tanai, = tan — 

t2+ti 2 



71 TT" 



(6.4) 



The hopping part of the Hamiltonian is now diagonal: 



(6.5) 



fc,cr 
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with 

e{k) = ^tl + tl + 2tit2 cos A; . (6.6) 
The Fermi velocity vp '^^ given by 

vp = , = , , (6.7) 



where Ad = 2{ti — t2) is the dimerization gap and W = 2{ti + is the total bandwidth 
(see Section. I). 

The Fourier transform of the Hamiltonian (|1.3|) describing the interaction is 



^ ki,k2,k3,k4 



where ki + k2 — — = Q is a vector of the reciprocal lattice. Usually Q = 0, but if one 
of the bands is half filled, the umklapp processes with Q = ±27r become important. Using 
the operators defined by Eq. (|6.3| ), the interaction can be written as 



Hint = -j^ L cosf + j)<AA,A4,T' (6-9) 

ki,k2,k3,k4 

where we have kept only the fermions near the Fermi energy. 

From the interaction part of the Hamiltonian we can identify the interactions between 
the electrons near the Fermi surface, the so called gf-couplings of the (^-ologyil. Usually these 
couplings are spin dependent. However, in our model the interaction is isotropic, thus we 
do not have to worry about the spin dependence, and they read: 

U 

91= 92 = 9i = ^ , 

■ (0 2 U 

,3 = -sm(2M- = -^YT(AWW 2 • ^'-'^^ 

The umklapp scattering amplitude 5^3 vanishes linearly with A/j for small A^/W , in agree- 
ment with the estimateJii obtained using perturbational arguments to get the strength of 
the umklapp scattering. Clearly, the model of Eq. ( |6.1(J| ) is not equivalent to the half-filled 
Hubbard model, for which all the g couplings are equal and given by U . 

14 



In perturbation theory, the logarithmic corrections to the vertex generate the differential 
equations of the renormalization group (RG) approach when one integrates out the degrees 
of freedom far from the Fermi leveiS. There are four differential equations altogether, but 
near half-filling we need only the equations which describe the charge degrees of freedom 
and give rise to the charge gap. As we know from Larkin and Sak0, to account correctly for 
the fluctuations which give rise to the ^JTj factor in Eq. (|5.4|) , we need the RG equations up 
to third order in (yf-s. Introducing g = 2g2 — Qi and denoting by g the renormalized couplings 
corresponding to momentum cutoff D, the 3^'^ order RG equations we need are 

do 1 2 1 ~~2 

^3 + 7^-i-i^^3 . 



d\nD TiVc' 27r2w2 



c 



dfi-s 1 ~- , 1 \~2~ , ~3l 

99i + -r^ 9 9^ + 9z , (6.11) 



d In D TTf c A'k'^vI 
where Vc = vp + (74/271 is the velocity of the charge excitations. There is a scaling invariant 
C = {g"^ — ^|)/(27rt>c — g) , and the system of equations above can be rewritten as a single 
differential equation: 

-^ + -^-— 1-7^:- • (6-12) 



dhiD V^^^c ^^"^c T^Vc) \ 2nvc 

Reducing the cutoff D, we scale towards the strong coupling region, where the RG 
equations are no longer valid. This cross-over occurs when the cutoff D corresponds to an 
energy scale that has been identified with the charge ga So integrating the scaling 

equations (|6.12|) , we get 



,nnD^^r ^ ( 4,, (l- . (6.13) 

where Dq is the initial cutoff with g{Do) = 9 a-nd g-^lDo) = g^. Keeping only the leading 
and next to leading terms of g, and introducing the notation C, = gs/g, one gets 

1 A ~ 1 /vri;, 1\ Arthy/T^ ^ 1 , ^ 

In lnDo = 7^1n 7 r, ^ +- + \nCA, 6.14 

where the constant IuCa contains the terms like l/g{Ac) which are small compared to 1/g, 
as g scales towards the strong coupling limit. As in that region the RG equations are not 
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any more valid, all the possible corrections are incorporated in that constantllll. We have 
separated 1/4 from the constant InCA for convenience. 

Replacing t>c by + 5'/27r { = g in the leading order in U), and exponentiating, we 

get 



Ann ^3 ( 7rt;pArthv^r^\ / lArthv^T^ l\ 

Its behavior in the case of small dimerization is given by 

A. = C.Do«./^eV'(f) . (6,16) 

where we have used the fact that Arth-\/l — ^ In (2/^) if ^ —>■ 0. This is the same form 
as found by LutherEl for the charge gap if the umklapp scattering is small compared to the 
2(72 — 91, which in our notation means ^ ^ 1 (see also Ref. [ 0). 

Using the g-couplings of Eq. (|6.10|) in Eq. ( |6.15|) provides us with the functional depen- 
dence of the gap on the repulsion U. This expression also includes a dependence on ^2/^1- It 
does noi give the full dependence on t2/^i however. The reason is the following: Our estimate 
of the charge gap accounts neither for the curvature of the band (as the RG equations are 
based on a model with linear dispersion relation) nor for the possible processes involving the 
empty band. This virtual processes will contribute as higher order corrections in U to the 
effective interactions near the Fermi surface. Since in the formula for the charge gap these 
effective couplings appear essentially as exp{-~TTVF/g), the U'^ correction to the effective cou- 
pling will give an 0(1) correction, i.e. if we assume for a moment that g = U — 'jU'^ + 0{U^) , 
then 1/g = 1/U + 7 + 0{U) and exp{—nvp/ g) = exp(— 7rt>j77) exp(— vrfj^/f/). So the cor- 
rection appears as a multiplying factor that depends on t2/ti in the equation for the gap. 
Clearly, to get a quantitative estimate of the gap one has to include these contributions. 

To be fully consistent with our determination of the U dependence, which relies on RG 
equations up to third order in g (see Eq. ( |6.11[ )), one should in principle calculate the vertex 



corrections to the coupling constants to third order in U. The problem with that program 
is that the calculation of the third order is hopelessly cumbersome. We think however that 
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a calculation of the vertex correction to second order in U is sufficient. The procedure we 
have used is the following: First, we calculate the second order corrections to the effective 
interactions near the Fermi surface by integrating out all the electron states except those 
which are closer to the Fermi surface than some small cutoff Dq in the lower band. Then, 
if the cut-off Dq is small enough, the dispersion relation of the electrons within this cut-off 
around the Fermi surface is essentially linear and we can safely use the RG equations of 
Eq. (|6.11|) to decrease the cutoff further down to Ac/vc- Now, the cut-off Dq can be taken 
very small: It just has to be larger than Ac/vc, and we already know that this quantity is 
exponentially small in the weak coupling limit. So this procedure should be valid. Finally, 
we will see that for small Dq (W/vc ^ Dq ^ Ac/vc) the result for the gap is independent 
of the cutoff Dq. 

To calculate the higher corrections in U, we must take into account the virtual processes 
that involve states in the upper band and consider the full interaction Hamiltonian instead 
of Eq. (^): 

jj Ul ^ f Ok, + - - ak^) , Q\ 

Hint = -J^ 2^ COS + - 

X (4i,T42a^fc3,i^fc4,T + fLjlAA,^ + dLA2jksJk,,^ + • • •) ' (6-17) 

where we have written only the terms which are responsible for the second order corrections 
in U. To that order, the effective interactions are given by 

U f/2 

g,{DQ) = - - — {D, + D,) , 
92{DQ) = ^-^iD, + D2 + Ds) , 

g,(DQ) = -smi2aF)j~^iD, + D,) , (6.18) 

where we have denoted hj Di, D^ the diagrams contributing to the effective interaction 
(see Fig. P). They are given by 



T7/2+D0 2tt 2eF — 2e{q) Jo 2% 2e{q) — 2eF 
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D2 



= 2 



= 4 



27r 
(ig 1 



e{q) - e{7r - q) 



27r2e^ + 2e(g) 



7r/2-Do _ cOs{aF - (ttg + Q;7r-5)/2) sin(Q;i7' + {aq + aT,-q)/2) 



271 



2eF + e{q) + £(7r - g) 



io 



dg — cos(q; — {aq + ctTr-g) / 2) sin(Q; J? + (dig + a^-q) /2) 



27r - £(7r - g) 

Here ap and s^;' stands for a{kF) and e{kF)- 



(6.19) 



The integrals and contains a logarithmic singularity ln(2/i)o), and, after 

some algebra, one can separate these contributions, so that the effective interactions can be 
written as 



^i(^o) 
^2(^0) 



U_ _ 1 

2 4 tivf 

U 1 f 1 - sin2(2Q;F) , 2 A 
^-^ — — In ^ + — 



4 ttvf 



-sin2(2a^)| [l+{A) - I, (A)] - ^ sm\2aF)IiiA) + Opo) 



gsiDo) = -sin(2Q;i7) 



U f/2 1 



4 TTVi? 



In^ - 

Da 



A 



I^iA)-I,~iA) 



^h{A) + P2iA) + OiDo) 



(6.20) 



where A — 2tit2/{t\ + and /i and I2 are non-singular integrals given by: 

r7r/2 



/■7r/2 1 

•^0 1 + a/i ± 



A (A) 



7r/2 



dtp 



+ y'l ±^cos((^) 

1 



y^l — Acos(</7) y^l + Acos((/?) / 

•^0 V Jl- 



y'l - C0S2((/?) / 2 + + A cos((/7) + ^1- A cos((/7) 



(6.21) 



For A ^ 1, they read: 



im = ^ T I + o(A^) , 
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h{A) = A + OiA'), 

/2(A) = 7r/4 + 0(^2) . (6.22) 

If it were only for a single band with linear dispersion relation these integrals near the 
ln(2/Z}o) would not appear. 

Using the couplings of Eq. (|6.20| ) in Eq. (|6.15|) , we get 



X exp (-i^ + i) e«<^)|l + 0{D,)] , (6.23) 

where, to first order, Dq disappeared from the expression of the gap. This expression now 
describes the full functional dependence of Ac on U/ti and t2/ii- The only thing that 
remains is to determine the prefactor Ca. This can be done as follows: From the Bethe 
ansatz solution of the Hubbard model, we know the value of the gap exactly in the limit of 
large dimerization [see Eq. (|5.4|) ]. It is easily checked that in the limit ^2 ^ the functional 
dependence of Eq. ( |6.23| ) is the same as in Eq. (|5.4| ). So we can use the result of Eq. ( ^.4| ) 



to determine the prefactor, and this gives Ca = '^y'^/'^- Our final expression for the gap in 
the weak coupling limit is thus 

/ lArthA l\ rtA) 
xexpf--— ^ + -je^(^), (6.24) 

where A = 2tit2/ {tl + 1^) and C{A) is given by 

C{A) = ^ [2(1 - A')I,{A) + 2A'I^{A) - (1 - 2A')h{A) - 2h{A) 
Arth /I r 

+ [-2Io{A) + (1 - A')h{A) + 2(1 - A')h{A)] . (6.25) 

It is interesting to compare this approach with that of Larkin and Sa k0 , who used slightly 
different arguments when they calculated the gap in the negative U Hubbard model in the 
small U limit. It can be shown very easily that our method is equivalent to theirs and gives 
the same result for the attractive Hubbard model. 
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B. Alternating on site energies 



The Fourier transform of the kinetic Hamiltonian ( |1.4|) using the Fourier transforms of 
the operators defined in Eq. ( |6.1|) is 

Ho = -eo E (4,.c,,. - - 2t cos(A:/2) ^ {e-''H,.Ck,. + ^''^.c,^.) ■ (6-26) 

fc,cr k,a 

It can be diagonahzed in terms of new creation and annihilation operators defined by 

dk,a = Ck,a COS Pk + ^''''^^5^^^ siu , 

f.a = c..sin/3, - e-^'^/^c.^cos/Jfc , (6.27) 



where jSk is given by 

tan 2Pk = - cos - , ^ ^ 4 j ' ^^'^^^ 

The kinetic part of the Hamiltonian reads 

i/o = - E <k) (4, A,. - fUk,.) , (6.29) 



with 



e{k) = ^Jel + 2t^ + 2t^cosk . (6.30) 
The dimerization gap opening at the Brillouin zone boundary is A^:) = 2eo, the 'total' 



bandwidth is W = 2Jel + At'^ and the Fermi velocity is 



vf = , . (6.31) 



Replacing the new operators defined in Eq. ( 6.27| ) for the lower and upper band into the 
interaction Hamiltonian ( |6.8|) we get: 

H:nt = U^ E [ 

{cospkr cos/3fc2 cospks cospki ± sin/5fcj sin/Jfc^ sinp^s sin/3fcj 4i,T42,i'^fe3,i'^fc4,T 
+ (sin cos (3k2 cos cos (3k4, =F cos (3ki sin /J^^ sin /?fc3 sin /3fcJ 
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+ (cos /3fc, cos /3fc2 sin /^fcg sin ± cos cos (3k2 sin /3fc3 sin (3k J 

^ {fL,lfL,ldk3,idk4,'\ + '^l4,T'^l3a^fc2,i.^fci,T) 
+ ...], (6.32) 

with the constraint ki + k2 — — = Q and the upper (lower) signs stand for the normal 
Q = (umklapp, Q = ±27r) scattering, respectively. There are altogether 16 terms, but we 
have kept only the terms which we need to calculate the g^s near the Fermi surface up to 
second order in U. 

Like in the case of alternating hopping amplitudes, to determine the effective couplings 
g, we first integrate out the high energy processes up to some cutoff Dq: 

giiDo) = [/(cosVf + sm^pp) - {Di + 2^2 + + 2 A) , 

^2(^0) = f/(cos^ Pf + sin^ Pf) - {Di + 20^ + D^, + D, + 2D^) , 

g^{b^) = f/(cos^ Pf - sin^ Pf) - {2Dj + Ds + Dg + 2D,o) , (6.33) 

where we have denoted hj Di, Diq the diagrams contributing to the effective interactions 
(see Fig. P) which are given by 

r'^ da ( cos^ Pf cos^ /3„ + sin^ Pf sin^ /3„ 

Di = 2 

J-k/2+Do 27r 2eF — 2e(g) 

r-n/2-bo dq (cos^ Pf cos^ Pq + siu^ Pf sin^ Pq^ 



2-K 2e{q)-2eF 

r'^ da ( cos^ Pf cos P„ sin /5„ — sin^ Pf sin /5„ cos /3„ 
D2 = 2 

J 1^/2 2tx 2eF 

■TT dq (cos^ Pf sin^ Pq + sin^ cos^ P^ 



D3 = 2 I ^ — 

Jo 2n 2eF + 2e{q) 

f'^/'^ dq (cos^ Pf sin Pq cos Pq — sin^ sin Pq cos 
^ ~^ '0 2^ 2l(^^ ' 

tt/2-Do dq 2 cos* Pf cos^ Pg cos^ /^Tf-q + 2 sin^ /Sf' sin^ Pq sin^ 
27r e(g)-£(7r-g) 
dq 2 CDS'* Pf cos^ Pg sin^ /Jth-^ + 2 sin'^ Pfsiio^ Pq cos^ /^jr-g 
27r e(g)+e(7r-g) 



D5 = -2 

1^6 = - 



'0 

■tt/2 
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D7 



dq cos^ Pf cos^ Pg sin^ /^jr-g — sin'* Pp sin^ Pg cos^ /3,r-g 



7r/2 27r 



2^^;' - + e{7i - q) 



dq cos"^ Pp sin^ /5q sin^ P-K-q — sin'* cos^ Pq cos^ /^jr-q 



27r 



2eF + £:(g) + e(vr - q) 



-4 



7r/2-Do cos"^ Pf cos^ /?g cos^ - sin^ /?F sin^ Pq sin^ 



27r 



^^72 dq cos^ Pf cos^ /J^ sin^ /5,r-g — sin** /^^ sin^ Pq cos^ 



(6.34) 



/o 27r + £(7r - q) 

Again, like for the case of alternating hoppings, we separate the logarithmic divergencies 
from integrals Di, and Dg, and we write the effective couplings as 



p2 I f2 Tj2 



g2{Do) 



el + 2t2 ttvf 

(I 



el + e\\ 2 ■ 

-4 ^ n ^ + Ji 

el + 2t^) Do 



ei + e 



U 



In ^ + /2 



£§ + 2^2 vrt;^ V 2 (£2 + 2^2)2 



£^0 



£^ + 2t2 



-.U - 



TTVf 



(el + 2t2)3/2 + 3^ 



(6.35) 



The expressions for the J's are rather complicated and, for brevity, we will give here only 
their value in the limit of large dimerization (the compounds which can be described by this 
model presumably have parameters which lie in this limit): 

32 + 137r f 



h 



32 



+ Oif/el) 



Ael 32 e^o 
t2 48 + Stt 

^2 d 



8 /^S\ 



44 



+ o(t74). 



(6.36) 



We have calculated numerically the Is for a few values of eo/t and they are presented in 
Table |IV|. 

For Qi, it is sufficient to take the value without the f/^ corrections: 



^ el + 2t2 



(6.37) 



Replacing the effective interactions into Eq. (|6.15|) , we get 



-t 



Ueo 



71 \2{6l + 2t2) 



exp 



^^el + 2t2Arth 



el + t\ 



exp 



1 t^ + el 
4 4^2 



Arth 



el + t^ 



(6.38) 
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where the constant has been fixed so that for large dimerization {eq ^ t) we get Eq. { f>.8\j . 
Here C is given by 

^2 



2^2 
o) 



1-^Arth ' 



~ (2^2 + £2)3/2 
-'3 12 



t2 



t4 



(6.39) 



VII. LARGE U LIMIT 

In the large f/ limit a canonical transformation can be applied to obtain an effective 
Hamiltonian analog to the t — J Hamiltonian of the non dimerized Hubbard model. We 
know already that if there are two alternating hoppings, a dimerization gap A/) opens at 
the Brillouin zone boundary, and, in the limit U +oo, this gap becomes the charge gap. 
This will be modified by the spin interaction, which will give corrections of order 1/U. 



A. Alternating hoppings 

If f/ ^ ^1,^2, and since we have two different hopping amplitudes, the effective model 
will be a t — J model with alternating ti and ^2 hoppings and alternating Ji and J2 exchange 
interactions: 

Htj = -ti E ^(4c,+i,, + h.c.)P-t2 E ^(4q+i,, + h.c.)P 

icvcn, cr iodd, cr 

+ Ji E {SiSi+i - ^niUi+i) + J2 E('S'i5'j+i - in^nj+i) 

icven iodd 
+ ^ E^(cU4+l,-.C,+i,,C,,+2,-. - cl,A+l,-^^^+l,-^^^+2,. + h-C)^ + ^(^7^^') • (7-1) 

Here the projector V = Yiil^ ~ ""-j,?^*,!] ensures that there are no doubly occupied sites, 
ni = nj^i + rii^i, and the generated exchange couplings are Ji = 4:t1/U, J2 = '^t\/U and 

To go further, let us follow the scheme developed by Shiba and Ogata in their study of 
the correlation functions for the large U Hubbard modefl. The basic idea is to treat the 
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exchange part as a perturbation. This is possible if the exchange integrals are smaller than 
the hopping integrals, which reduces to the condition Ji <^ t2 (remember that t2 < h, and 
thus that J, J2 < Ji), or, in terms of the original parameters, U ^ 4t^/t2- In this limit, 
the part describing essentially hopping of spinless fermions can be solved exactly and the 
exchange part is then treated as a perturbation. This can be achieved by assuming that 
the wave function is the product of a charge and a spin wave function, where the spin wave 
function is defined in a Hilbert space of dimension 2^, i.e every charge has the additional 
freedom of having its spin up or down: 

1^) = ® 1$) . (7.2) 

Here |$) is the spin part of the wave function, while the wave function of the N spinless 
fermions, is the ground state of the kinetic Hamiltonian. Following standard perturbation 
technique for the case of a degenerate ground state, we will have to diagonalize the following 
2^ X 2^ Hamiltonian to hft the degeneracy of |$) and get the 1/U energy corrections: 

(HuY/L^ -t,{ciciy-t2{c\c2y 

(J J \ 1 ^ 1 

-j{non,y + -fin.n^y - J(4clc,c,)'j - E(5,5,+i - -) , (7.3) 

where the (A)' denotes the expectation value of the operator A in the Fermi sea of spinless 
fermions, i.e. (^4)' = {sf\A\sf). 

Before we continue with the calculation of the gap, we need to know the expectation 
values of electron operators in the Fermi see of spinless fermions, where kF,sf — ^kp- Using 
the translational symmetry of the system, it is enough to calculate the following basic 
expectation values 

{clc,y = jY.^''ncWy. 
^ k 

^ k 

(4c,r=yEe^'^'/'(40', (7.4) 
^ k 
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where j is even everywhere. More comphcated expressions can be calculated by using Wick's 
theorem, i.e. by constructing all possible pairings. For example: 



{clc\c,c,y = n{clc,y - {clc,y{c{c,y . (7.5) 

In the following, the sums above will be replaced by integrals, 

2 1 



The free fermion Hamiltonian ( |6.2D has been solved in Sec. VI. Using the diagonal electron 
operators d and / defined in Eq. (|6.3|), we get 



(cicfc)' = ^^^ind,k - rif^k) , (7.7) 

where n^ fc = ci^cifc, similarly for ray ^, and the spin index of the fermion operators has been 
dropped. For j even, this gives 



smiTjn 
nj 



(7.8) 



This is the same expression as for a non-dimerized model. For {c[cj)', we get: 

1 f-^vo 



27r Jo 

where is defined in Eq. (|6.4|) and 



{c\c,y = — / rfA;cos(A;(j - l)/2 - a,) , (7.9) 



Tin , if 72 < 1/2: 

^o={ (7.10) 
7r(l -n) , if n > 1/2. 



After some algebra, we get 



TT JO 



/ t \/ 1 tiCOS(^J +t2COS(^ J - 2 

(cic,) = z I dip — ; ^ ^ • (7.11) 

tl + tl + 2tit2Cos2ip 
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For j = and j = 2, this can be written 

1 



{clcoY 

{c\c,y 



2nti 
1 

2^ 



[{h-t2)F{^o,q) + {ti + t2)E{^o,q)] 
[{t2 - ti)F{ipo, q) + (ti + t2)E{^^, q)] 



(7.12) 



where we have introduced the notation 



q 



ti + h 

and where the eUiptic integrals E{ipo, q) and -^(v^o, q) are defined by 



(7.13) 



^(<^o,g) =- I — J 

JO V 1 "~ q sm yj^ 
E{'fo,q)=J^ dv^y^l - g^sinv?^ 



(7.14) 



Now, let us turn back to the eff'ective Hamiltonian (|7.3|) . For the spin degrees of freedom, 
we are left with an A^-site Heisenberg Hamiltonian with an effective exchange interaction 
Jeff which we can read from Eq. ( [7.3|) : 



JeS = y(no'^l)'+ y(^1^2)'- -/(CqCiCiCs)' 



2 
U 



2!^2 , ^ ^ sin 27m (^1-^2) ^2/ ^ 

^ -f^ wo^q) 



n^(ti +^2) -ntits 



(7.15) 



where we have used Eqs. ( |7.5| ) and (|7.12| ) 

In the ground state of the AF Heisenberg model the expectation value (^KSoS*! — |)|$) 
is — In 2 in zero magnetic field, as shown by Griffithsill. The energy per site (\E'|if(j|\E')/L is 
then 



£tj — £sf — Jeff In 2 



(7.16) 



where £sf = ~'ti{ciCo)' — t2{c\c2)' is the kinetic energy of the spinless fermions per site. From 
Eq. ( fml) , we get 

tl + t2 



TT 



'-E{(po,q) 



(7.17) 



which in the non-dimerized case [q = 1) gives the correct energy (2t/7r) sinvm. 
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For n = 1/2, we get 



tf + t2 {ti — t2 



-K\q) 



(7.18) 



where K{q) = F(7r/2,g) is the complete elhptic integral defined in Eq. (|7. 14|) . 

The energy Etj/L has a cusp as a function of filling at = 1/2, and this nonanaliticity 
comes from the (po [see Eq. ( [r.lOp ]. Therefore the first derivative of Etj/L has a jump at 
quarter filling, which is nothing but the charge gap: 



de 



tj 



dn 



de 



tj 



n=l/2+0 



dn 



(7.19) 



n=l/2-0 



In our case this gives 

Ae = 2{ti - t2) 



^ 41n2ti + t2.. /2v^ 
1 — — A 



TT U 

where we have used the following identities: 



ti + h 



(7.20) 



dn 
dE{ipo,q) 



dF{^o,q) 



dn 



ri.=l/2+0 



n=l/2+0 



dn 
dE{ipo,q) 



2tt 



dn 



n=l/2-0 



n.=l/2-0 



2tiJI - g2 



(7.21) 



Using the limiting behavior of the elliptic function K{q) = ln(4/ y/1 — g^), the gap in the 
small dimerization limit is 



' 41n2ti+t2, 4(ti + t2) 
1 In 

TT U tl-t2 



(7.22) 



or 



A,. = A 



D 



2\n2W , m' 
1 m 



(7.23) 



7r U Ad_ 

In Fig. |1^ we have compared the analytical results presented above with the numerical 
data. The agreement is already very good at U/ti = 10 for ^2/^1 between 0.3 and 0.5. To 
get such a good agreemnt for ^2/^1 = 0.1, one has to go to larger values of U, so that the 
condition U ^ 4:tl/t2 ~ 40 is satisfied. For large values of t2, the gap is small and the 
numerical estimate is not accurate. 
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As t2 0, Eq. (|7.20|) gives Ac = 2ti — (41n2)tf/f/, in apparent contradiction with 
the second of equations (|4.1|), which gives Ac = 2ti — Atl/U for t2 = 0. The origin of 
the discrepancy is that Eq. ( |7.20| ) holds if f/ ^ ^1/^2, while the second of equations ( |4.1| ) 
requires that ^2 ^ tf/U- These conditions are clearly incompatible. Actually the cross-over 
between these two regimes can be observed on Fig. |T0|, where the points below the solid line 
{t\/U ~ 0.1) gives Ac ~ 2ti — 4t1/U, while for larger values of U we get the correct 4 In 2 
slope. 



B. Alternating on site energies 

Like in the previous case, a canonical transformation can be applied in the limit U ^ t, eg 
to obtain a t — J like effective Hamiltonian 

i^cr 2evcn,(T 

i 

+ ^ E^(cU4+i,-.c,+i,.c,+2,„, - clA+i,-.c,+,,-.'C,^2,a + h.c.)P + 0(1' /u') , (7.24) 

where the projector V = Yiil^ — '>^i,i'>^i,i\ ensures that there are no doubly occupied sites and 
the generated exchange coupling is J = At^ /U. If the additional condition J ^ cq is satisfied 
(that is, if f/ ^ 4t^/eo), we can again follow Ogata and Shiba. 

Let us now calculate the expectation values for this model. The free fermion model is 
given by Eq. ( |6.26|) and it has been solved using transformations Eq. (|6.27|) . So, for the 



expectation values we get 



(45.)' = ^^^^^^ - ^^^^^^ 
(4c.)' = e--/^^^^^^^^sin2A 

(4c,/ = e-^/^ ^^'^~^^'^ sin2/j,. (7.25) 
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Using the results we obtained in Eq. ( |7.8|) , for j even we get 



J- , smn jn 



sin TTjn 



TTJ 



(7.26) 



and after some algebra 



(4c_i)' = (4c)' = — (Jj + /j_2) 
^0 

where the Ij denote the following integral: 

^ 2 27r 72^.0 



(7.27) 



1 r'i'Po 



27r Jo 



dk 



cos{kj/2) 



^l + {2t/eo)^cos^{k/2) 



(7.28) 



and ifio is defined by Eq. (|7.1CI|) . These are elliptic integrals and e.g. Iq and I2 are given by 



^0 



^(<^o,g) 



where 



2 ^oaAo + 

- /2 = -^^2 [^(^0, q) - E{^o, q)] , 



2t 



(7.29) 



(7.30) 

The energy of the system is Etj = (\E'|iJjj|\l') and the energy per site etj = Etj/L is still 
given by Eq. ( [7.16| ), but now the kinetic energy of the spinless fermions per site is 

sM = 2t(4c,)' - 1 ((4co)' - {c\c,y) 

E{^o,q), (7.31) 



el + 4t2 



where we have used Eqs. (|7.27| ) and ( [7.29| ). For the effective exchange coupling we get 



Jeff — J 



{noriiY - - {{clc\c^C2y + {c\clc2C^y) 



2 

sin 2iTn 



— <n— n 

U 1 27r 



2tiH^ 



(7.32) 
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where we have used Wick's theorem, namely: 



(C1C2C2C3)' = ^o(cIc3)' - (cjci)'^ 



and 



(?^o?^i)' - 2 ((C0C1C1C2)' + (cIclcaCg)' 



noui - ^(rai(4c2)' + rao(clc3)') 

2 sin 2?™ 
n -n— loilo-h) . 



For n = 1/2, we get 

Je^ = ^\t'- ^K{q) [K{q) - E{n/2, q)] | 
The gap, using Eq. ( |7.19| ), is 



Ae = 2eo{l 



21n2 V^o + 4t2 
~ U 



[l + q')K{q)-E{n/2,q) 



and in the small dimerization limit it is given by 
A, = 2eo 



1 — — m — 



TT 



(7.33) 



(7.34) 



(7.35) 



(7.36) 



(7.37) 



U eo 

If we parametrize the gap in this limit with and W, we can see that it is identical with 



Eq. {^^. 
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FIGURES 

FIG. 1. Schematic representation of the model with: a) alternating hopping amplitudes ti 

and t2; b) alternating on-site energies with energy splitting 2eo . 

FIG. 2. Band structure of the model with U = 0. 

FIG. 3. Schematic picture of the different regions as a function of model parameters for the 
model with: a) alternating hopping amplitudes; b) alternating on-site energies. 

FIG. 4. a) Approximation for the charge gap for t2/ti =0.1 to 0.9 in increments of 0.1 from the 
top to the bottom. The diamonds are the numerical data for ^2/^1 =0.1 to 0.5; b) Approximation 
for the gap, but now only from t2Ai =0.6 to 0.9 in increments of 0.1 from the top to the bottom. 

FIG. 5. Numerical estimates of the gap (diamonds) for t2 = 0.1 compared with the analytical 
results. 

FIG. 6. Weak coupling expression compared to numerical estimates for t2 = 0.3 (upper curve) 
and t2 = 0.4 (lower curve). 

FIG. 7. Numerical values of the gap Ac(L) for L=4,8,12 and 16 as a function of 1/L for 
U/ti = 5. Note that Ac{L) converges very well for small values of t2/ti. 

FIG. 8. Perturbational corrections to the effective backward scattering gi, forward scattering 
g2 and umklapp scattering 53 for the model with alternating hoppings shown up to the second 
order. The solid (dashed) line represents the fermions in the lower d (upper /) band. 

FIG. 9. Perturbational corrections to the effective backward scattering gi, forward scattering 
g2 and umklapp scattering g^ for the model with alternating hoppings shown up to the second 
order. The solid (dashed) line represents the fermions in the lower d (upper /) band. 



33 



FIG. 10. Numerical estimates of the gap in the large U limit (diamonds) together with the 
analytical approximation. The results presented are for t2 =0.1, 0.2, 0.3, 0.4, 0.5 and 0.7 from the 

top to bottom 
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TABLES 

TABLE I. The parameters of the gap equations for the weak coupUng and strong couphng 
hmit for the model with alternating hoppings. 



t2/tl 


a 


b 




c 


0.1 


0.36887 


0.6336 


1.8 


1.6816 


0.2 


0.46782 


1.2990 


1.6 


2.0167 


0.3 


0.49990 


2.0304 


1.4 


2.3984 


0.4 


0.48604 


2.8669 


1.2 


2.8368 


0.5 


0.43672 


3.8588 


1.0 


3.3474 


0.6 


0.36100 


5.0790 


0.8 


3.9555 


0.7 


0.26867 


6.6518 


0.6 


4.7075 


0.8 


0.17046 


8.8399 


0.4 


5.7054 


0.9 


0.07709 


12.4449 


0.2 


7.2658 




TABLE IL 


The parameters Uq, ai, 


bi and 62. 





t2/tl 


Uo/h 


ai 




b2 


0.2 


1.30 


0.45312 


1.5438 


13.0786 


0.3 


2.03 


-0.13375 


1.6173 


13.5738 


0.4 


2.87 


-0.97220 


1.5540 


10.7001 


0.5 


3.86 


-1.91825 


1.4287 


3.3094 


0.6 


4 


-1.84743 


2.6357 


0.0736 


0.7 


4 


-1.50580 


5.3353 


-2.2707 


0.8 


4 


-1.09045 


11.5364 


-7.4996 


0.9 


4 


-0.60222 


33.3189 


-34.0969 
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TABLE III. The parameters of the gap equations for the weak coupHng and strong couphng 

limit for llic iiKKicl willi allernaliug on site cuergieti. 



eo/t 


a 


b Ao/t 


c 


0.2 


0.60292 


8.82123 0.4 


5.63087 


0.4 


0.78553 


6.00854 0.8 


4.41422 


0.6 


0.90155 


4.53739 1.2 


3.71014 


0.8 


0.97686 


3.61670 1.6 


3.21882 


1 


1.02292 


2.98899 2 


2.84615 


2 


1.03186 


1.56009 4 


1.78537 


3 


0.93780 


1.04544 6 


1.28117 


4 


0.84788 


0.78494 8 


0.99167 


5 


0.77446 


0.62816 10 


0.80616 


TABLE IV. The Is for some values of eo/t. 


eo/t 


h 




h 





(1/8) In 2 








0.2 


0.0832216 


0.0007691 


0.0020951 


0.4 


0.0739705 


0.0026394 


0.0039468 


0.6 


0.0613605 


0.0046787 


0.0053306 


0.8 


0.0480877 


0.0061374 


0.0061143 


1 


0.0360938 


0.0067624 


0.0063123 


2 


0.0066400 


0.0038914 


0.0035338 


3 


0.0013135 


0.0015298 


0.0014380 


4 


0.0003309 


0.0006442 


0.0006189 


5 


0.0001034 


0.0003044 


0.0002961 
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